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Abstract. We first review the known mathematical results concerning the 
KP type equations. Then we perform numerical simulations to analyze various 
qualitative properties of the equations : blow-up versus long time behavior, 
stability and instability of solitary waves. 



1. Introduction 

The main goal of this paper is an analysis of qualitative properties of solutions to 
generalized Kadomtsev-Petviashvili (KP) equations via a numerical approach. 

The (classical) Kadomtsev-Petviashvili (KP) equations 

(1) {Ut + UxxX + UUj;)x ±Uyy =0 

were introduced in [57] to study the transverse stability of the solitary wave solution 
of the Korteweg- de Vries equation which reads in the context of water-waves 

(2) Ut^Ux^ uux + (^ - T)uxxx = 0, X G M, t > 0. 

Here T > is the Bond numbeij^ which measures surface tension effects in the 
context of surface hydro dynamical waves. 

Actually the (formal) analysis in [57] consists in looking for a weakly transverse 
perturbation of the one-dimensional transport equation 

(3) Ut^Ux = 0. 

This perturbation, which is obtained by a Taylor expansion of the dispersion rela- 
tion uj{ki,k2) = ^/kl -\- kl of the two-dimensional linear wave equation assuming 
<C 1 and <C 1, amounts to adding a nonlocal term, leading to 

(4) Ut^Ux^ ^d~^Uyy = 0. 

Here the operator is defined via Fourier transform. 



Date: October 28, 2010. 

''^Note that here T = corresponds to the absence of surface tension. In the Fluid Mechanics 
community, the Bond number is often defined as the inverse of our T. 
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= ^/(O, where ^ = 
Remark 1. Equation Q is reminiscent of the linear diffractive pulse equation 

2utx 

where Ay is the Laplace operator in the transverse variable studied in [8]. 

The same formal procedure is applied to the KdV equation assuming that the 
transverse dispersive effects are of the same order as the x-dispersive and nonlinear 
terms, yielding the KP equation in the form 

(5) Ut^Ux^ uux + - T)uxxx + ^d~^Uyy = 0. 

By change of frame and scaling, ([5| reduces to ([T]) with the + sign (KP II) when 
T < I and the - sign (KP I) when T > |. 

Note however that T > ^ corresponds to a layer of fluid of depth smaller than 0.46 
cm, and in this situation viscous effects due to the boundary layer at the bottom 
cannot be ignored. One could then say that "the KP I equation does not exist 
in the context of water waves" , but it appears naturally in other contexts (see for 
instance Remark 2 below). 



Of course the same formal procedure could also be applied to any one-dimensional 
weakly nonlinear dispersive equation of the form 



(6) Ut^Ux^ f{u)x - Lux = 0, X G M, t > 0, 

where f(u) is a smooth real- valued function (most of the time polynomial) and L 
a linear operator taking into account the dispersion and defined in Fourier variable 

by 



(7) Lum=pmu{£,), 

where the symbol p(^) is real-valued. The KdV equation corresponds for instance 
to f{u) = and p(^) = — Examples with a fifth order dispersion in x are 
considered in [6], j^, [59] . 

This leads to a class of generalized KP equations 



(8) Ut-^Ux^ f{u)x - Lux + -d^ ^Uyy = 0, X G M, t > 0. 



Let us notice, at this point, that models alternative to KdV- type equations ([6| are 
the equations of Benjamin-Bona-Mahony (BBM) type [10] 

(9) Ut^Ux^ uux + Lut = 

with corresponding two-dimensional "KP-BBM-type models" (in the case p{^) > 0) 

(10) Ut -\- Ux -\- uux + Lut + d~^dyU = 
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{Ut + 



'X 



+ uux + Lut)x + dlu = 



S{t) 



It was only after the seminal paper [57 that Kadomtsev-Petviashvili type equations 
have been derived as asymptotic models to describe the propagation of long, quasi- 
unidirectional waves of small amplitude (under an appropriate scaling) in various 
physical situations. This was done for instance formally in [3 in the context of 
water waves (see [74], [75] for a rigorous approach in the same context). 

A rigorous approach (with error estimates) was also made by Ben-Youssef and 
Lannes [llj for a class of hyperbolic quasilinear systems. We refer also to [ST, [9^ 
for a rigorous derivation from a Boussinesq system or the Benney-Luke equation. 
In all these works, the convergence rate is proven to be low (we will be a bit more 
precise in Subsection 2.1), a phenomena which is deeply analysed in [74 (see also 
[75] ) and which is mainly due to the singularity at = of the KP dispersion 
relation. 

On the other hand, Dryuma [26] found a Lax pair to the KP- I/II equations, 
proving the "integrability" of the KP equations (see [97J for a precise description 
of the Inverse Scattering aspects of the KP equations). 

Remark 2. In some physical contexts (not in water waves!) one could consider 
higher dimensional transverse perturbations, which amounts to replacing d~^Uyy 
in (|8| by 9~^A^ix, where is the Laplace operator in the transverse variables. 

For instance, the KP I equation (in both two and three dimensions) also describes 
after a suitable scaling the long wave transonic limit of the Gross-Pitaevskii equation 
(see [13 for rigorous results for the solitary waves (ground states) in 2D and [23] 
for the Cauchy problem in dimensions two and more). 

More precisely, let the Gross-Pitaevskii equation be 



The Gross-Pitaevskii equation can be written in hydrodynamic form via the Madelung 
transform provided does not vanish. This is the case in the transonic regime where 
in particular is close to one. In this regime, the KP I equation describes, after 
a suitable rescaling, the behavior of = 1 — 

Note again that in the classical KP equations, the distinction between KP I and 
KP II arises from the sign of the dispersive term in x, which characterizes what 
physicists call "positive" or "negative" dispersion media. 

A noteworthy fact implied by the derivation of KP equations is that, as far as the 
transverse stability of the KdV solitary wave '^1-soliton^' 



(12) 



itpt + AV^ + (1 - |V^P)V^ = 0, in x R, 



for functions ip with finite Ginzburg- Landau energy 



(13) 




(14) 
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is concerned, the natural initial condition associated to ([T]) should be of the type 
'^0 = + ^0 where vq is either "localized" in x and or localized in x and 
^/-periodic. 

The same observation is of course valid for (|8| in connection with the transverse 
stability of solitary waves of (|6|. 

This paper is organized as follows. In the next Section we review the main rigorous 
mathematical results concerning KP type equations. They are mainly obtained by 
PDE techniques with the exception of results obtained via the Inverse Scattering 
machinery for the classical KP I and KP II equations. This Section has an interest 
in itself but it leads also to conjectures on the long time dynamics (or blow-up in 
finite time when it is expected) for KP type equations. To keep this survey short 
we will only discuss the main results and ideas and refer to the original papers 
for details and proofs. Also we will mainly comment on results obtained by PDE 
or Nonlinear Analysis techniques. We refer for instance to the excellent survey 
article [27 which describes in particular results on the Cauchy problem obtained 
by Inverse Scattering methods for the classical KP I and KP II equations. 

The aim of the next Section 3 which is the heart of this work is to give numerical 
evidence of those conjectures and to suggest further theoretical investigations. We 
will also provide precise numerical decay rates of various norms of the solutions. 

Notations 

For any real number sq, the notation Sq means "for any s > sq". The norm in the 
classical Lebesgue spaces L^(R^) will be denoted by \.\p. 

We denote, for any real number 5, by i^^(R^) the Sobolev space of distributions 
u in the Schwartz space ot tempered distributions <S'(M'^) equipped with the norm 
||i^||5 = |(l + |^P)^'u|2 where u denotes the Fourier transform of u. 

We will use sometimes the anisotropic Sobolev spaces H^^'^'^(M?)^ equipped with 
the norm 

IHU,.. = l(i + l6n*(i + l6P)*«|2. 

For the homogeneous version of those spaces, we just replace the weight (1 + 
ICP)* by 



2. A SURVEY OF THEORETICAL RESULTS ON KP EQUATIONS 

We survey here various mathematical results concerning KP type equations. They 
will serve as a guide for our numerical simulations, together with the (many) open 
problems. 

2.1. The zero- mass constraint in x. In (|8|, it is implicitly assumed that the 
operator d~^dy is well defined, which a priori imposes a constraint on the solution 
which, in some sense, has to be an x-derivative. This is achieved, for instance, if 
u G 5'(R2) is such that 

(15) ^i'e2u{t,^u^2)eS\R^), 

thus in particular if ^ S(t, ^i, ^2) ^ <S'(]R^). Another possibility to fulfill the 
constraint is to write u as 

d 

(16) u{t,x,y) = --v{t,x,y), 
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where v is a continuous function having a classical derivative with respect to 
which, for any fixed y and t ^ 0, vanishes when x ±00. Thus one has 



(17) 



u{t^ X, y)dx = 0, G M, t 7^ 0, 



in the sense of generalized Riemann integrals. Of course the differentiated version 
of ([8|, namely 

{Ut ^Ux^ UUx - LUx)x + dyU = 0, 



(18) 



can make sense without any constraint of type (15) or (17) on and so does the 
Duhamel integral representation of (|8|, 

(19) u{t) = S{t)uo- [ S{t- s){u{s)ux{s))ds, 



where S{t) denotes the (unitary in all Sobolev spaces i^^(R^)) group associated 
with (|8]), 

(20) S{t) = 6-*^^--^^-+^-"'^') . 

In particular, the results established on the Cauchy problem for KP type equations 
which use the Duhamel (integral ) formulation (see for instance |22] |lQ8j ) are valid 
without any constraint on the initial data. 

One has to be careful however in which sense the differentiated equation is taken. 
For instance let us consider the integral equation 



(21) u{x,y,t)=S{t)uo{x,y)- [ S(t - t')[u{x,y,t')ux{x,y,t)]dt' , 

Jo 

where S{t) is here the KP II group, for initial data uq in i^^(R2), 5 > 2, (thus uo does 
not satisfy any zero mass constraint), yielding a local solution u G C([0, T]; H^(M?)). 



By differentiating (21 ) first with respect to x and then with respect to t, one obtains 
the equation 



dtdxu ^ dx{uux) ^ d^u ^ d^u = in C([0,T];i7^ 



However, the identity dfdxU = dxdtu holds only in a very weak sense, for example 
in D'((0,T) X M?). 

On the other hand, a constraint has to be imposed when using the Hamiltonian 



formulation of the equation. In fact, the Hamiltonian for (18) is 

^2 



(22) 



and the Hamiltonian associated with (11) is 



(23) 



It has been established in [Q^ that, for a rather general class of KP or KP-BBM 



equations, the solution of the Cauchy problem obtained for (18), ( |TT| ) (in an appro- 
priate functional setting) satisfies the zero-mass constraint in x for any t 7^ (in a 
sense to be precised below), even if the initial data does not. This is a manifestation 
of the infinite speed of propagation inherent to KP equations. Moreover, KP type 
equations display a striking smoothing effect (different from the one reviewed in 
Subsection 2.3 though) : if the initial data belongs to the space L^(R^) ni^^'°(R^) 
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and ifixGC([0,T];i72'0(R2))[2] 

is a solution in the sense of distributions, then, for 
any t > 0, u{.^t) becomes a continuous function of x and y (with zero mean in x). 
Note that the space L^(R^) fi i^^'^(R^) is not included in the space of continuous 
functions. 

The key point when proving those results is a careful analysis of the fundamental 
solution of KP-type equations which turns out to be a x derivative of a continuous 
function of x and with respect to x, which, for fixed t ^ and y^ tends to 

zero as X ^ ±00. Thus its generalized Riemann integral in x vanishes for all values 
of the transverse variable y and of t 7^ 0. A similar property can be established for 
the solution of the nonlinear problem [93 . Those results have been checked in the 
numerical simulations of [69 , [68j. 

We refer [30], |112j for a rigorous approach of the Cauchy problem with initial data 
which do not satisfy the zero-mass condition via the Inverse Spectral Method in 
the integrable case. 

Nevertheless, the singularity at ^1 = of the dispersion relation of KP type equa- 
tions make them rather had asymptotic models. First the singularity at ^1 = 
yields a very bad approximation of the dispersion relation of the original system 
(for instance the water wave system) by that of the KP equation. 

Another drawback is the poor error estimate between the KP solution and the 
solution of the original problem. This has been established clearly in the context 
of water waves as we briefiy recall now. 

Let us first recall the KP scaling for water waves. We denote by a a typical ampli- 
tude, h the mean depth, A (resp. ji) the typical wavelength in x, (resp. y). Then 

we set 

where e <Cl, 5'i - 1, ^'2 - 10 

It has been established rigorously in [74], [75] that the KP II equation yields a 
poor error estimate when used as an asymptotic model of the water wave system. 
Roughly speaking, the error estimates with the relevant solution of the (Euler) 
water waves system reads: 

\\UEuler - UkpW = o{l) . 

So the error is o(l), {0{^/e) with some additional constraint) instead of O(e^t), 
which should be the optimal rate in this regime (as it is the case for the KdV, 
Boussinesq equations or systems). Nevertheless the KP II equation reproduces 
(qualitatively) observed features of the water wave theory. For instance the well- 
known picture below displays the interaction of two oblique "line solitary waves" 
on the Oregon coast which shows a striking resemblance with the so-called KP II 
2-soliton. 

We refer to [75] for the derivation of a new class of weakly transverse approximations 
of the water waves system in the Boussinesq regime which are purely local (without 
the anti-derivative in x), have a dispersion relation which fitting well with that of 



We will see in Subsection 2.4 that KP type equations (in particular the classical KP I and 
KP II equations) do possess solutions in this class. 

^ In the case of KP II, one can use the explicit form of the fundamental solution found in lOOJ. 
e is thus the small parameter which measures the weak nonlinear and long wave effects. 
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Figure 1 . interaction of line solitons. Oregon coast 

the original Euler system, and which provide optimal error estimates. No mass 
constraint has to be imposed there. 

2.2. Localized solitary waves. Sohtary waves are solutions of KP equations of 
the form 

where (= y or (y^z)) is the transverse variable and c > is the solitary wave 
velocity. 

The solitary wave is said to be localized if tpc tends to zero at infinity in all directions. 
For such solitary waves, the following energy space is natural 

(with an obvious modification in the three dimensional case), 

and throughout this section we will deal only with finite energy solitary waves. 

Due to its integrability properties, the KP I equation possesses a localized, finite 
energy, explicit solitary wave, the lump [84^ : 



(24) 



(pcix - ct,y) 



8c(l - f(a;-ct)2 + |-y2) 
(1 + f(a;-ct)2 + ^fy2)2 • 



Another interesting explicit solitary wave of the KP I equation which is localized in 
X and periodic in y has been found by Zaitsev [128J . It reads 



(25) 
where 



Zc{x,y) = 12a^ 



1 — /3cosh(ax) cos{6y) 
[cosh(ax) — (3 cos{Sy)]'^ 



{a,l3) G (0,(X)) X (-1,+1), 



8 



C. KLEIN AND J.-C. SAUT 



and the propagation speed is given by 

24-/3^ 

Let us observe that the transform a ia^ 5 ^ iS^ c ^ ic produces solutions of the 
KP I equation which are periodic in x and locahzed in y. 

We address first the question of non existence of solitary waves for the following 
classes of KP type equations 

(26) ut + u^Ux + Uxxx + ^d~^Uyy = 0, 
where e = ±1, in the two-dimensional case, and 

(27) Ut + u^Ux + Uxxx + (id~^Uyy + hd~^Uzz = 0, 
where a, 6 = ±1, in the three-dimensional case. 

By establishing suitable Pohojaev type identities, it has been proved in [T9^ that 
(under mild regularity conditions), no localized solitary waves exist when c < in 
both dimensions and for c > when 

e = —1, p > 4 or e = 1 and p arbitrary, 

ab = —1 (resp. a = b = 1) and p arbitrary, 
or 



a = b = -1, p> -. 

In particular, the 2D KP II- type equations {e = 1, p > 1) do not possess any 
localized solitary waves. This was precised by de Bouard and Martel [18j who 
proved that the classical KP II equation does not possess L^- compact solutions 
(those finite energy solutions are more general than solitary waves). 

We refer to [19] for nonexistence results of localized solitary waves of KP type 
equations with a cubic-quintic dispersion in x. 

We turn now to the existence of localized, finite energy solitary waves. We will 
consider as a basic example a class of generalized KP I equations in two dimensions 
of the form 

(28) Ut + u^Ux + Uxxx - d~^Uyy = 0, 

where p = ^, n odd, m and n relatively prime. 

Solitary waves are looked for in the energy space Y'(R^) which can also be defined 
(see [19]) as the closure of the space of x derivatives of smooth and compactly 
supported functions in M? for the norm 
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The equation of a solitary wave of speed c is given by 

(29) cd,i^ - ^d^^ - dl^ + d-\dl^) = 0, 
which imphes 

(30) cd,,^ - {^d,^), - dt^ + dl^ = 0, 



When e F(]R^), the function d^^dyil) is weh-defined (see [19]), so that (29) makes 



y 

sense. 



Given any c > 0, a sohtary wave ?/^c of speed c is deduced from a soUtary wave 
of velocity 1 by the scaling 

(31) i^c{x,y) = c^lJl{^/cx,cy). 



We now introduce the important notion of ground state solitary waves. 
We set 



^ JR2 ^ Jr2 Z[P^ Z) J ^2 



2 ^ 7r2 ^KP^ ^) 7r2 

and we define the action 



S{N)=Ekp{N)^^ [ N\ 
^ Jr2 

We term ground state^ a solitary wave N which minimizes the action 5* among all 
finite energy non-constant solitary waves of speed c (see [19 for more details). 

It was proven in [19 that ground states exist if and only if c > and 1 < p < 4. 
Moreover, when 1 < p < |, the ground states are minimizers of the Hamiltonian 
Ekp with prescribed mass (L^ norm). 

Remark 3. When p = 1 (the classical KP I equation), it is unknown (but conjec- 
tured) whether the lump solution is a ground state. 

In the three dimensional case, ground states exist when a = b = —1 and 1 < p < | 
[19]. They are never minimizers of the corresponding energy with fixed norm 
(actually this minimization problem has no solutions in this case). 

It turns out that qualitative properties of solitary waves can be established for a 
large class of KP type equations. Ground state solutions are shown in [20 to be 
cylindrically symmetric^ that is radial with respect to the transverse variable up to 
a translation of the origin. 

On the other hand, any finite energy solitary wave is infinitely smooth (at least 
when the exponent p is an integer) and decay with an algebraic rate in two 
dimension and r~^+^ in dimension three [20 . Actually the decay rate is sharp in 
the sense that a solitary wave cannot decay faster that r"'^, d = 2, 3. 

Moreover a precise asymptotic expansion of the solitary waves has been obtained 
by Gravejat ([3j]). 
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The fact whether the ground states are or are not minimizers of the Hamiltonian 
is strongly linked to the orbital stability of the set Qc of ground states of velocity 
c. Note that the uniqueness^ up to the obvious symmetries, of the ground state of 
velocity c is a difficult open problem, even for the classical KP I equation. 

Saying that Qc is orbitally stable in Y means that for 0c ^ Qc^ then for all e > 0, 
there exists S > such that if uq G Y]^ is such that ||^io — 0c||y ^ then the 
solution u{t) of the Cauchy problem initiating from uq satisfies 



sup inf \\u{t) — ^p\\Y ^ ^^ 
t>o i^^O 

Actually it is proved in pT that the ground state solitary waves are orbitally 
stable in dimension two if and only if 1 < p < | while there are always unstable in 
dimension three. 

We will see below the link between stability and blow-up in finite time of the 
solutions to the Cauchy problem. 

Remark 4. No rigorous stability (or instability) result seems to be known for the 
Zaitsev solitary wave solution of the KP I equation. 



2.3. The linear group. As previously mentioned, the linear part associated to 
any KP type equations (18) generates a group S{t) = ex.p{t{dxL + d~^dyy)) which 
is unitary in all Sobolev spaces, in particular in L^. For instance in the case of KP 
I/II equations, we have 



Sxpi/iiit) = exp{t{da:xx T dyy)) 

When the problem is set in R^, Sxpi/iiit) has nice dispersive properties. In fact 
an analysis of the fundamental solution leads to the dispersive estimate |lQ3j 

C 

(32) ||^kp////(^)</>||l-(r2) < ^||(/>||li(r2), Vt ^ 0. 

Such an estimate leads to L^L^y Strichartz estimates [103" which are useful to solve 
the nonlinear Cauchy problem. Typically, one gets estimates of the type (here S{t) 
stands for the KP I or the KP II group) : 



(33) ||5'(^)0||l'2(R;L-(r2) < C||0||l2(r), 

for any r G [2, +oo), and where q = q{r) is defined by ^ = | 
Similar estimates hold for the "Duhamel integral" 



A/(t) 



fs{t 

Jo 



s)f{s)ds. 



Of course an estimate like ( 32 ) is false when the spatial domain is the torus or 
the "mixed " case M x T or T x R. 



On the other hand, dispersive estimates of the type ( [32| hold true when the Uxxx 
dispersive term in the KP equation is replaced by dxP{dx)u^ where P is an even 



'Some extra regularity on uq is actually needed. 
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polynomial (see [9^), in particular for the "cubic-quintic" KP equations. In this 
case the corresponding Strichartz estimates involve a global smoothing effect [9 . 

Another interesting dispersive property of the KP group is the local smoothing 
property. 

It has been known since the seminal paper of Kato [60 that the Airy group, that 
is S{t) = e*^^^, satisfies the following local smoothing estimate 

l|5x5'W^)||L2(0,T;Lf^^(R) < C'l | (/)| | (R)) • 

In other words there is (locally in space) a gain of one derivative. 

A similar property was established in |lQ3j for the KP II group. Actually, for any 
initial data (j) G L^(R^), d~^(l) E I/^(M^), the corresponding solution u of the linear 
KP II equation satisfies the estimate, for any T > and > 0, 



J-T J-R J-oo 

which also displays a local gain of one derivative. 

A stronger smoothing effect holds for the nonlinear KP II equation when the initial 
data decay fast to the right [79^. More precisely, if the initial data satisfies 



for all positive integers L, then the corresponding solution u (which exists globally, 
see the next Section) satisfies u{.^t) G C^(]R^) for all positive t's. 

A result of this type (with different hypothesis on the initial data) was established 
in [78] for the nonlinear KP I equation. 

2.4. The nonlinear Cauchy problem. All the KP type equations can be viewed 
as a linear skew-adjoint perturbation of the Burgers equation. Using this structure, 
it is not difficult (for instance by a compactness method) to prove that the Cauchy 
problem is locally well-posed for data in the Sobolev spaces H^^ 5 > 2 in both 
dimensions (see |124j , [103] , [50] for results in this direction). 

Unfortunately, this kind of result does not suffice to obtain the ^/o6a/ well-posedness 
of the Cauchy problem. This would need to use the conservation laws of the equa- 
tions. For general KP type equations, there are only two of them, the conservation 
of the norm and the conservation of the energy (Hamiltonian). For the general 
equation ([8| where f{u) = ^^i^^^^, and without the transport term Ux (which can 
be eliminated by a change of variable), the Hamiltonian reads 




JR2 



< +00 



(34) 



uLu + {d^ ^'^yY + 



p + 2 



and for the classical KP I/II equations 



(35) 




where the + sign corresponds to KP I and the — sign to KP 11. 
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Note that the "integrable" KP I and KP II equations possess more conservation 
laws, but not infinitely many as it is often claimed (see below). 

In any case, it is clear that for KP II type equation, the Hamiltonian is useless to 
control any Sobolev norm, and to obtain the global well-posedness of the Cauchy 
problem one should consider solutions, a very difficult task. On the other hand, 
for KP I type equations, one may hope (for a not too strong nonlinearity) to have 
a global control in the energy space, that is 

(36) E = {ue uLu e d'^Uy e L^}. 

For the usual KP I equation, E = Y^ which was defined in the previous subsection. 

The problem is thus reduced to proving the local well-posedness of the Cauchy 
problem in spaces of very low regularity, a difficult matter which has attracted a 
lot of efforts in the recent years. 

Remark 5. By a standard compactness method, one obtains easily the existence of 
global weak finite energy solutions (without uniqueness) to the KP I equation (see 

At this point, there is a striking difference between the KP I and the KP II equation, 
linked to their different dispersion relations. The KP II equation can be solved by a 
Picard scheme implemented on the integral (Duhamel) formulation, implying that 
the flow map uq u{.^t) is smooth in natural Sobolev spaces. In this sense the KP 
II equation can be viewed as a semi-linear equation. 

On the other hand, it has been proved in [105 for the periodic case and in [92] 
for the Cauchy problem in M? that the KP I equation cannot be solved by such 
a procedure, in any natural Sobolev class. This implies that the flow map cannot 
be smooth, and actually it cannot be even uniformly continuous on bounded sets 
of the natural energy space Y (see [71 ). This is a typical property of quasi-linear 
hyperbolic systems {eg the Burgers equation). It is thus natural to consider the KP 
I equation as a quasi-linear equation. 

Those ill/well posedness by iteration properties have a huge importance when solv- 
ing the Cauchy problem since the methods are quite different in both cases. We 
will first focus on the classical KP I and KP II equations. 

For the KP II equation, a breakthrough was made by Bourgain f22^ who proved 
that the Cauchy problem for the KP II equation is locally (thus globally in virtue 
of the conservation of the norm) for data in L^(]R^), and even in L^(T^). This 
result was later improved ([m], [HT], [122] [118], [53], [39], [40]) to allow the case 
of initial data in negative order Sobolev spaces. 

Remark 6. Bourgain's proof uses in a crucial way the fact (both in the periodic 
and full space case) that the dispersion relation of the KP II equation induces the 
triviality of a certain resonant set. A similar property was used by Zakharov [130J 
to construct a Birkhoff formal form for the periodic KP II equation with small initial 
data. On the other hand, the fact that for the KP I equation the corresponding 
resonant set is non trivial is crucial in the construction of the counter-examples of 
|105j . [92] and is apparently an obstruction to the Zakharov construction for the 
periodic KP I equation. 

The KP II equation in three space dimensions has been studied in |121j . [52] . The 
Cauchy problem is locally well-posed for initial data in anisotropic Sobolev spaces 
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having 1+ derivative in x and 0+ derivatives in z in L^(R^). It is an open problem 
whether or not those solutions are globally defined. 

We refer to [103], [50], [39], [65], [66], [35], [4^ for similar results on generalized 
KP II type equations. Note that for nonlinearities higher than quadratic, that is 
f{u) = p > 1 in ([8|, and Lu = —u^^^ one gets only the local well-posedness. 

Neither global existence nor blow-up in finite time are established. 

On the other hand, global results are known to be true for higher order KP II 
equations, eg for the fifth order equation corresponding to Lu = —Uxxxx [106 . 

Most of those results are based on iterative methods on the Duhamel formulation, 
using the Fourier restriction X*'^ spaces of Bourgain, or variant of them, possibly 
with the injection of various linear dispersive estimates (see jl03j . [9]). 

We also would like to mention the paper by Kenig and Martel [6T where the Miura 
transform is used to prove the global well-posedness of a modified KP II equation. 

In view of [92 the situation is quite different for the KP I equation. Note however 
that the KP I equation with a fifth order dispersion in x is globally well-posed in 
the corresponding energy space [107 (see also [3T for local well-posedness results 
without the "zero mass constraint" on the initial data). 

For the classical KP I equation, the first global well-posedness result for arbitrary 
large initial data in a suitable Sobolev type space was obtained by Molinet, Saut 
and Tzvetkov [94 . The solution is uniformly bounded in time and space. 

The proof is based on a rather sophisticated compactness method and uses the 
first invariants of the KP I equation to get global in time bounds. It is worth 
noticing that, while the recursion formula in [131 gives formally a infinite number of 
invariants, except for the first ones, those invariants do not make sense for functions 
belonging to L^(R^) based Sobolev spaces. 

For instance, the invariant which should control | \uxxx{'-, ^)| |l2 contains the L^ norm 
of d~^dy{u^) which does not make sense for a non zero function u in the Sobolev 
space H^{M?). Actually (see [93], [94]), one checks easily that if d~^dy{u'^) G 
L2(R2), then J^^dy{u'^)dx = dy j^^u^dx = 0, Wy e M, which, with u e ^^(R^), 
implies that u = 0. Similar obstructions occur for the higher order "invariants". 

One is thus led to introduce a quasi- invariant (by skipping the non defined terms) 
which eventually will provide the desired bound. There are also serious technical 
difficulties to justify rigorously the conservation of the meaningful invariants along 
the ffow and to control the remainder terms 

The result of [94 was extended by Kenig [62^ (who considered initial data in a 
larger space), and by lonescu, Kenig and Tataru [49] who proved that the KP I 
equation is globally well-posed in the energy space Y. 

Remark 7. For both the KP I and KP II equations, the Inverse Scattering method 
provides only global well-posedness for small initial data. We refer to \127\ for the 
KP II equation and to ^132^, [112^ for the KP I equation. 

Concerning the generalized KP I equations 



(37) ut + Uxxx - ^Uyy + u^Ux = 0, 

an important remark is that the anisotropic Sobolev embedding (see [T2^) 
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/ iwr2dxd2/<ciiziiij''iiu,ii^.p-iw,ii|„ 

which is vahd for < p < 4, imphes that the energy norm \\u\\y is controhed in 
term of the norm and of the Hamiltonian 



(p+l)(p + 2) 



dxdy, 



if and only if p < | . 

This suggest global well-posedness when p < | and blow-up when p > |, a fact 
which will be proved in the next Section. 

2.5. Blow-up issues. We consider here a class of generalized KP equations 



(38) ut + Uxxx + ^Uyy + u^Ux = 0, 

where e = 1 for the KP II equations and —1 for the KP I equations. 

The KP I equations are focusing while the KP II are defocusing and to this as- 
pect they share some properties with the focusing (resp. defocusing) nonlinear 
Schrodinger equation (NLS). In particular one could expect the blow-up in finite 
time of the solutions of generalized KP I equations with a sufficiently high nonlin- 
earity. 

For the focusing NLS, one can prove the blow-up in finite time of the norm of 
the spatial gradient by using a virial identity (see eg |lllj ). For the generalized 
KP I equation, one can use a "transverse" virial identity to prove the blow up 
in finite time of the norm of the transverse gradient (that is to say in 2D of 
I !%(•? ^)| |l2(r2)). This identity has been formally derived in pr20j and rigorously 
justified in |lQ3j (one needs to prove a local well-posedness theory in a weighted 
space). It reads for solutions of the generalized KP equations in a suitable 
functional class: 

(39) ^I{t)=4e[-pE{t)^^ [ uidxdy ^ ^^^e f {d'^Uyfdxdy], 

(40) Fit) = Jj-iul - eid-^u,f - ^—^^]d.dy, 
and 

I(t) = / y^u^dxdy. 

Jr2 

A first direct consequence of (|39| for KP I equations (e = —1) is that, when the 
initial data uq (is "regular" and) belongs to the space 

S = {(/) G ^, # G L\R^), (/) / 0, E{^) < 0}, 

the corresponding solution of the Cauchy problem cannot remain smooth for all 
times if p > 4. More precisely, there should exist T* > such that 
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(41) ]hnJ\uy{.,t)\\L2 =+oo. 

A similar result is obtained in the three-dimensional case ie when d~^Uyy in (38) is 



replaced by ^A^u, where is the transverse Laplacian in the y and z variables. 



The critical value is then p = 2 and (41) should be replaced by 



(42) \\mJ\V^u{.,t)\\L2 =+oo. 

Remark 8. It is readily seen that the set S is not empty. It contains in particular 
functions of the type ±A9^exp(— — y"^) when A is a large positive number, the 
sign depending on the "parity" of p. 

Those results are not fully satisfactory since we already observed that the natural 
energy norm is controlled by he norm and the Hamiltonian if and only if < p < 
I in the two-dimensional case. In the three dimensional case, the critical exponent 
is p = I which suggests that the "standard" {p = 1) 3D- KP I equation is not 
supposed to have global solutions for arbitrary initial data. 

This result was improved by Liu [81 who used invariant sets of the generalized KP 
I flow together with the virial argument above to prove the existence of initial data 
leading to blow-up in finite time of \\uy{.^t)\\L2(j^2^ when p > |. He also proved a 
strong instability result of the solitary waves when 2 < p < 4. 

Similar results are established in [82] for the three-dimensional generalized KP I 
equation, in the range 1 < p < |. 

Remark 9. A detailed analysis of the blowing-up solution (profile, blow-up rate,..) 
is lacking. 

Remark 10. No rigorous finite time blow-up has been established for the periodic 
generalized KP I equation. 

Remark 11. No result seems to be known concerning the long time behavior (finite 
time blow-up versus global existence) of the solutions to the generalized KP II 



equations (38) when e = 1 and p >2). 



Remark 12. The finite time blow- up does not a priori invalidate the applicability of 
KP equations as asymptotic models. As all asymptotic models, they are supposed 
to describe the dynamics of a more general system, for a suitable scaling regime 
(here weakly nonlinear, long wave with a wavelength anisotropy) and on relevant 
(finite!) time scales. For instance, in the water wave context, the KP equations 
cease to be relevant models for times larger than O(^). 

2.6. The Cauchy problem in the background of a non localized solution. 

As was previously noticed, it is quite natural in view of transverse stability issues 
to consider the Cauchy problem for KP equations on the background of a solitary 
wave of the underlying KdV equation. 

The Inverse Scattering method has been used formally in [4] and rigorously in [5] to 
study the Cauchy problem for the KP II equation with non decaying data along a 
line, that is 2x(0, x, y) = Uoo{x — vy) + 0(x, y) with 0(x, 7/) ^ as + ^ oo and 
'^oo(^) ^ as |x| ^ oo. Typically, u^q is the profile of a traveling wave solution 
U (k.x — ujt) with its peak localized on the moving line k.x = ujt. It is a particular 
case of the N- soliton of the KP II equation discovered by Satsuma jlQ2j (see the 
derivation and the explicit form when N = 1,2 in the Appendix of [97 ). As for 
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all results obtained for KP equations by using the Inverse Scattering method, the 
initial perturbation of the non-decaying solution is supposed to be small enough in 
a weighted space (see [5j Theorem 13). A similar result has been established by 
Fokas and Pogrebkov [29] for the KP I equation. 

On the other hand, PDE techniques allow to consider arbitrary large perturbations 
of a class of non-decaying solutions of the KP I/II equations. 

We will therefore study the initial value problem for the KP I and KP II equations 

(43) {Ut + Uxxx + UUx)x ± Uyy = 0, 

where u = u{t,x,y)^ {x^y) G M^, t G M, with initial data 

(44) u{0, X, y) = (t){x, y) + V^c(^, v), 

where ?/^c is the profil^of a non- localized (i.e., not decaying in all spatial directions) 
traveling wave of the KP I/II equations moving with speed c 7^ 0. 

We recall that, contrary to the KP I equation, the KP II equation does not possess 
any localized in both directions traveling wave solution. 

The background solution ?/^c could be for instance the line soliton (1-soliton) of 



the Korteweg- de Vries (KdV) equation (14), or the N-soliton solution of the KdV 



equation. A/" > 2, or in the case of the KP I equation, the Zaitsev soliton. 

The KdV N-soliton is of course considered as a two dimensional (constant in y) 
object. 

This problem can be attacked in two different ways. Either one considers a localized 
(in X and y) perturbation, or the perturbation is localized in x and periodic in y. 

Both frameworks have been applied to the KP I and KP II equations, and they lead 
to global well-posedness of the Cauchy problem for arbitrary large initial perturba- 
tions. More precisely, for the KP I equation, the Cauchy problem with initial data 
which are localized perturbations of the KdV N-soliton or of the Zaitsev soliton is 
globallly well-posed [95 . The corresponding result for the KP II equation is proved 
in [961. 



The global well-posedness of the KP II equation with arbitrary initial data in (R x 
T) is proved in [96j. The corresponding result for the KP I equation with initial 
data in a suitable Sobolev space is established in [48] , 



2.7. Transverse stability issues. The KP I and KP II equations behave quite 
differently with respect to the transverse stability of the KdV 1-soliton. 

Zakharov [129 has proven, by exhibiting an explicit perturbation using the inte- 
grability, that the KdV 1-soliton is nonlinearly unstable for the KP I flow. Rousset 
and Tzvetkov |lQlj have given an alternative proof of this result, which does not 
use the integrability, and which can thus be implemented on other problems {eg for 
nonlinear Schrodinger equations). 

The nature of this instability is not known (rigorously) and it will be investigated 
in our numerical simulations. 

On the other hand, Mizomachi and Tzvetkov [89 have recently proved the L^(RxT) 
orbital stability of the KdV 1-soliton for the KP II flow. The perturbation is 
thus localized in x and periodic in y. The proof involves in particular the Miura 



*This means that ip(x — ct, y) solves |43 
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transform used in [64 to established the global well-posedness for a modified KP 
II equation. 

Such a result is not known (but expected) for a perturbation which is localized in 
X and y. 

2.8. Comparison between the KP and KP/BBM type equations. As pre- 
viously noticed, the KP/BBM equations ([lO| are relevant counterparts to KP type 
equations, for instance to the KP II equation 

The Cauchy problems for those regularized KP equations has been studied in [108 , 
for initial data in the natural energy space. We refer also to for more regular 
initial data and for a stability analysis of the solitary waves in the focusing case. 

The comparison of the dynamics of the KdV and BBM equation has been investi- 
gated both theoretically and numerically in [17 . A similar study was made in [86] 
between the KP and KP/BBM equations. Recah that KdV and BBM (resp. KP 
and KP/BBM) equations are supposed to describe the same dynamics, on relevant 
time scales. 

2.9. Long time dynamics. KP type equations can be classified into two cate- 
gories, the defocusing ones (KP II type), and the focusing ones (KP I type). For 
the former, the long time dynamics is expected to be governed by dispersion and 
scattering. In particular, at least for small initial data the sup norm should decay 
with the linear rate, that is 1/t (see |lQ3j ). 

On the other hand, the dynamics of the latter case is expected to be governed by 
the solitary waves and blow-up in finite time is expected in the supercritical case. 
To this extent the situation is reminiscent of that of the nonlinear Schrodinger 
equations |lllj . 

Relatively few mathematical results are know concerning those issues, and one 
goal of the present paper is to present numerical simulations in support of (or in 
motivation of) various conjectures on the long time dynamics. 

For the generalized KP equation 

Ut + Uj^XX ± dx^Uyy + U^Uj: = 0, 

it was proven in [43] that if p > 2 and if the initial data is sufficiently small, then 

ikwiioc < c(i + \t\)-\\og{2 + \t\r, \\u,m\oo < c{i + \t\)-\ 

for all t G M, where = 1, if p = 2 and = if p > 3. 
Moreover, a large time asymptotic for Ux can be obtained. 

Remark 13. No such result seems to be known for the KP II equation, p = 1 (see 
however Theorem 9.3 in [130 for a scattering result with small initial data by Inverse 
Scattering methods). We also mention [40] where a scattering result for small data 
in the scale invariant non isotropic homogeneous Sobolev space H~'^^'^'^{M?) is 
obtained. 

The large time behavior of solutions to the KP II equation with arbitrary large 
data is unknown. An interesting open problem is that of possible uniform in time 

^Note however that when the BBM trick is apphed to the KdV or to the KP equation with 
strong surface tension (Bond number greater than |) one gets a ill-posed equation, so that what is 
called KP-I/BBM equation is merely a mathematical object without, as far as we know, modelling 
relevance... 
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bounds on the higher Sobolev norms. It has been proven in [54j that the eventual 
growth in time of Sobolev norms is at most polynomial (see also |123j for long time 
bounds for the periodic KP II equation). 

The long time dynamics of KP I type equations for arbitrarily large data is also an 
important open problem. For the standard KP I equation, it is expected that the 
ground state solutions (whether or not the lump is one of them) should play a role 
in this dynamics as the KdV solitary waves for the corresponding equation. 

2.10. Initial data in the Schwartz class. One can prove that the Schwartz 
space <S(]R^) of functions with rapid decay is not preserved by the (linear and 
nonlinear) flows of KP type equations, in fact the rapid decay is not preserved as 
it is easily seen on the solution of the linear KP I or KP II equations which reads 
in Fourier variables 

^(6,6,^) = ^^o(6,6)exp(it(^? ± f )) 

Taking for instance i^o equal to a gaussian function, one sees that u{.^t) is not even 
in L^(R^) for t 7^ since this would imply by Riemann-Lebesgue theorem that 
u{.^t) is a continuous function of ^ which obviously is not the case. This proves in 
particular that the solution at time t cannot decay faster than ^ at infinity, where 
r = \/ + The same conclusion can be drawn in the nonlinear case. In fact 
the solution of the nonlinear equation with the same gaussian initial data has the 
Duhamel representation in Fourier variables 



^(ei,6,,t) = iio(ei,6)exp(zt(e?±f )) + ^ Aiexp(z(t-5)(^?±f ))^:2(^i,6,^)^^ 

K\ ^ Jo a 

and one checks that the integral part defines, for any fixed t > 0, a continuous 
function of ^ while the free part is not continuous as shown previously. 

A stronger decay of solutions is obtained when extra (rather unphysical) conditions 
are imposed on the initial data uq^ namely the vanishing of 'So(^i, ^2) at higher and 
higher order as ^1 ^0. 

The proofs above extend to any KP type equation like (pi), the crux of the matter 
being the singularity at ^1 = of the kernel 

The fact that the solutions of KP type equations do not decay fast at infinity yields 
difficulties discussed in Section 3.1 for the realistic numerical simulations of the 
solutions in the whole space. 

2.11. Varia. We mention briefiy here a few other mathematical questions concern- 
ing KP type equations. 

Many papers have been recently devoted to the study of initial boundary-value 
problems for various dispersive equations such as the KdV equation, motivated in 
particular by applications to control theory. 

A IBV problem for the KP II equation on a strip or a half space has been studied 
in [87]. We do not know of any similar result for the KP I equation. 

Many works have been recently devoted to the control or stabilization of linear or 
nonlinear dispersive equation. 
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We are not aware of any work on the control of (linear or nonlinear) KP type 
equations. A unique continuation property (which could be useful in control theory) 
for the KP II equation is established in [98 . A corresponding result for the KP I 
equation does not seem to be known. 

As was previously alluded to in the context of water waves, one sometimes need for 
physical reasons to add a viscous term to the KP equations leading to interesting 
mathematical problems. For a viscous term of the Burgers type, that IS '^XX 7 
Molinet and Ribaud [91 have established the global well-posedness of the Cauchy 
problem for data in L^(R^) for the KP-II-Burgers equation and in the energy space 
Y for the KP-I-Burgers equation. This result has been extended by Kojok [77 in 
the case of the KP-II-Burgers equation for initial data in the anisotropic Sobolev 
space H^^^{E?), for any s > -\. 

A sharp decay asymptotic analysis of the solutions of KP-Burgers equations as 
t +00 is carried out in [90]. 



3. Numerical study of generalized KP equations 

In this section we will numerically study the stability of exact solutions to the KP 
equations for various perturbations, and the blow-up phenomenon in generalized 
KP equations. 



3.1. General setting. We will now study numerically the solution to initial value 



problems for the generalized KP equations (38). We will only consider functions 
u that are periodic in x and y with period Lx27r and L^27r respectively or in the 
Schwartzian space of rapidly decreasing functions, which are essentially periodic for 
numerical purposes if restricted to an interval of length L on which the functions 
decrease below values that can be numerically distinguished from zero (10~^^ in 



double precision). As noted in subsection 2.1, all solutions of (38 ) satisfy in this case 
for t > the constraint /^^^ Uyydx = 0. To obtain solutions which are smooth at 
time t = (for convenience we only consider times t > 0) , the initial data also have 
to satisfy this constraint. It is possible to numerically solve initial value problems 
which do not satisfy the constraint, see for instance [68 , but the non-regularity in 
time is numerically difficult to resolve. Therefore we will only consider initial data 
that satisfy this condition, typically by imposing that uo{x^y) is the x-derivative 
of a periodic or Schwartzian function. This even holds for perturbations of exact 
solutions we will consider in the following. 

The imposed periodicity on the solutions allows the use of Fourier spectral methods. 



Equation (38) is equivalent to 



(45) dtU^^^^-if^u^X^^^u = 0. 

The singular multiplier is regularized in standard way as — i/(^i + ^0). Nu- 

merically this is achieved by computing — i/(^i + ie) where e is some small number 
of the order of the rounding error. Since we only consider functions u satisfying the 
constraint ^2^(t,^i = 0,^2) = 0, this does not lead to any numerical problems. 



For a numerical solution equations (45) are approximated by a truncated Fourier 
series, i.e., a discrete Fourier transform for which fast algorithms exist. With this 
approximation one obtains a system of ordinary differential equations for which an 
initial value problem has to be solved. This system is stiff due to the high order 
spatial derivatives, especially since we want to resolve strong gradients. Standard 
explicit schemes would require prohibitively small time steps for stability reasons in 
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this case, whereas unconditionahy stable imphcit schemes are typicahy not efficient, 
especially if one aims at high precision. In [67] we compared several fourth order 
schemes for the KdV equation, and in [TCT a similar analysis was performed for KP. 
It was found that exponential time differencing schemes are most efficient in this 
case, see for instance the reviews [88], [44]. There are several such schemes known 
in the literature which show similar performance. We use here the fourth order 
scheme proposed by Cox and Matthews [25 . 

As noted in subsection |2.1Q[ solutions to generalized KP equations for Schwartzian 
initial uq in general will not stay in the Schwartz class, unless uq satisfies in addition 
an infinite number of constraints. For generic initial data, the solutions u{t) will 
develop 'tails' decaying only algebraically in certain spatial directions. This implies 
for the used periodic setting that Gibbs phenomena will appear on the borders of 
the computational domain for larger times. These jump discontinuities affect the 
numerical accuracy. Thus a much higher resolution than for KdV (where solutions 
to Schwartzian initial data stay in the Schwartz space) would be needed to obtain 
high numerical precision. In the below computations this effect is addressed by 
choosing a larger computational domain than would be necessary for similar KdV 
computations. The algebraic tails will lead due to the periodicity condition to 
'echoes' which are, however, in most cases of comparatively small amplitude. 

The numerical accuracy of the computations is controlled by propagating exact solu- 
tions and by computing conserved quantities of KP which are numerically functions 
of time due to unavoidable numerical errors, see [70]. As discussed before gener- 
alized KP equations conserve the L^-norm (or mass) and the energy (40). Since 
the latter quantity contains the term d~^dyu(t^x^y) which implies in Fourier space 
a division by ^i, an operation that affects accuracy in the determination of the 
quantity, we will trace mass conservation. In [70 it was shown that relative mass 
conservation A := |1 — M(t)/M(0)| (M being the numerically L^-norm) overesti- 
mates typically the accuracy of the numerical solution in an sense by 2 orders 
of magnitude. Thus we always aim at a A smaller than 10~^ to ensure accuracy of 
plotting precision. Obviously the quantity A only makes sense if the spatial numeri- 
cal resolution is of the same order, i.e., if the Fourier coefficients S(t, ^i, ^2) decrease 
for large ^1,^2 to the same order of magnitude. The precision of the computations 
will therefore always be controlled by mass conservation and Fourier coefficients. 
In addition we check whether the found numerical results remain the same within 
the considered accuracy if computed with higher or lower resolution. 

Exact solutions to the KdV equation are automatically solutions to both KP equa- 
tions without explicit ^/-dependence. Thus KdV solitons are line solitons of KP, 
i.e., solutions that are exponentially localized in one spatial direction and infinitely 
extended in the other. As mentioned above the KdV soliton (14) is unstable for 
KP I and lineary stable for KP II. To test the code we first propagate initial data 
from the KdV soliton uq = Usoi{x — xq) for c = 4 and xq = —2Lx {L^ = 8 and 
Ly = 8^ X e [— TT, 7r]I/a^, y G [—7r^7r]Ly) with = 2^ and Ny = 2^^ Fourier modes 
with Nt = 3200 time steps from t = until t = 6. For both KP I and KP II we 
find at t = 6 that the relative mass conservation A = 3.89 * 10~^ and that the L"^- 
norm of the difference between numerical and exact solution is equal to 9.96* 10~^. 
This shows that the code is well able to propagate the exact solution even for the 
unstable case. It also illustrates that the quantity A can be used as expected as an 
indicator for the numerical accuracy. In the following subsections we will test how 
periodic perturbations of order unity of the KdV soliton are propagated by the KP 
equations. 
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3.2. Perturbations of the KdV soliton for KP II. The similarity between KP 
I and KP II for the exact KdV sohton initial data disappears if we perturb these 
data for instance by 

(46) 

Up = 6{x - xi) exp(-(x - xi)^) (exp(-(?/ + L^7r/2)^) + ex.p{-{y - L^7r/2)^)) . 

The perturbations are Schwartzian in both variables and satisfy the constraint ( |2.1[ ). 
They are of the same order of magnitude as the KdV soliton, i.e., of order 0(1), 
and thus test the nonlinear stability of the KdV soliton. 

Notice that the computations are carried out on a doubly periodic setting, i.e., 
and not on R^. Thus the perturbations cannot be radiated to infinity, but will 
stay in the computational domain. Since they are of a similar size as the initial 
soliton, they will lead to a perturbed soliton though the KdV soliton is known to be 
stable for KP II. One possibility is that the initial perturbation is essentially homo- 
geneously distributed over the whole computational domain. Another possibility 
would be that they form eventually a KdV soliton which is, however, not observed 
in the below examples. 

For KP II we use the same parameters for the computation as for the tests for the 
KdV soliton, but now for initial data uo{x^y) = 12sech^{x-\-2Lx) —Up {xi = xo/2), 
i.e., a superposition of the KdV soliton and the not aligned perturbation. At time 
t = 6 we obtain for the numerical conservation of the mass A = 4.29 * 10~^, almost 
the same precision as for the propagation of the soliton. The result can be seen in 
Fig. [2] The perturbation is dispersed in the form of tails to infinity which reenter 
the computational domain because of the imposed periodicity. The soliton appears 
to be unaffected by the perturbation which eventually seems to be smeared out in 
the background of the soliton. 

The situation changes slightly if the perturbation and the initial soliton are cen- 
tered around the same x- value initially, i.e., the same situation as above with 
Xi = xq = —2Lx' In Fig. [3] we show the difference between the numerical so- 
lution and the KdV soliton Usoi for several times for this case. It can be seen that 
the initially localized perturbations spread in ^/-direction, i.e., orthogonally to the 
direction of propagation and take finally themselves the shape of a line soliton. 
The perturbations appear to modulate the soliton here. It is difficult numerically 
to determine the long time behavior. Extending the computation till time t = 16, 
one finds a similar behavior as for t = 6 as can be seen in Fig. [4] Thus it appears 
that the perturbations of the KdV soliton travel here partly with the soliton and 
lead to oscillations around it. Since the speed of the soliton is not affected, it is 
still the original, not a soliton with slightly increased c which would change both 
height and speed of the traveling wave. But these perturbations traveling with the 
soliton seem to be present for very long times. 

A similar behavior is observed if the soliton shape itself is deformed, for instance 
in the form of the initial data uo{x) = 12sech^(x + 0Acos{2y/Ly)). In Fig. [sjthe 
difference of the solution for KP II for this initial data and the line soliton can be 
seen. The computation is carried out with A^^^ = 2^^, Ny = 2^ ^ = 16, Ly = S 
and Nt = 6400 with a final A - 10"'^. It can be seen that the solution approaches 
the KdV soliton with perturbations oscillating around it. This is even more visible 
in Fig. [6] where the L^-norm of the solution is given. It can be seen that the latter 
oscillates around the value for the KdV soliton. Numerically it is difficult to decide 
whether these oscillations will finally disappear. 

3.3. Perturbations of the KdV soliton for KP I. Perturbations of the KdV 
soliton are propagated in a completely different way by the KP I equation. It 
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t=2.1 




Figure 2. Difference of the solution to the KP II equation 
for initial data given by the KdV soliton plus perturbation, 
u{0^x^y) = Usoiix + 2La;,0) and perturbation Up = 6{x — 
xi) exp(-(x - (exp(-(y + Ly7r/2)^) + exp(-(^ - Ly7r/2)^)) , 
xi = —Lx and the KdV soliton for various values of t. 



is known that the line soliton is unstable for this equation. One of the reasons 
for this is the existence of lump solitons (24). Note that the lump soliton only 
decays algebraically in both spatial directions, as |?/| ^ oo. Because of this and 
the ensuing Gibbs phenomena in the periodic setting used here for the numerical 
simulations, the appearance of lumps will always lead to a considerable drop off in 
accuracy. In contrast to the line solitons it is localized in both spatial directions. 
Explicit solutions for multi-lumps are known, see [84 . It was shown by Ablowitz and 
Fokas [2 that solutions to the KP I equation with small norm will asymptotically 
decay into lumps and radiation. 

Exactly this behavior can be seen in Fig. [7| for the perturbed initial data of a line 
soliton with the perturbation (46) and xi = xq = —2Lx^ i.e., the same setting as 
studied in Fig. [3] for KP II. Here the initial perturbations develop into 2 lumps 
which are traveling with higher speed than the line soliton. The formation of these 
lumps essentially destroys the line soliton which leads to the formation of further 
lumps. It appears plausible that for sufficiently long times one would only be 
able to observe lumps and small perturbations which will be radiated to infinity 
if studied on M?. Notice that the algebraic fall off towards infinity causes the 
above mentioned numerical problems. This is reflected by the fact that the Fourier 
coefficients decrease much slower than in the case of KP II such that N^, = 2^ and 
Ny = 2^^ Fourier modes are needed to ensure that the Fourier coefficients go down 
by at least 5 orders of magnitude in and ^2- For the relative mass conservation 
one obtains A = 4.3 * 10"^ with N, = 6400. 
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Figure 3. Difference of the solution to the KP II 
equation for initial data given by the KdV soliton 
Usoi{x + 2Lx^0) plus perturbation Up = 6{x — xi)exp( — (x — 

xi)2)(exp(-(y + L^V2)')+exp(-(y-L^V2)')), = -2L„ 
and the KdV soliton for various values of t. 




Figure 4. Difference of the solution to the KP II equation for the 
situation of Fig. [3] and the initial KdV soliton for t = 16. 
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t=0.0 t=2.56 




Figure 5. Difference of the solution to the KP II equation for 
the initial data uo{x) = 12sech^(x + 0Acos{2y/Ly)) and the KdV 
soliton for several values of t. 

2.14 I ^ ^ ^ ^ ^ ^ ^ 1 
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Figure 6. L^-norm of the solution to the KP II equation for the 
situation of Fig. [5] in dependence of t. 
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t=0.0 



t=1.4 





t=2.8 



t=4.0 





Figure 7. Solution to the KP I equation for initial data given 
by the KdV soliton Usoi{x + '^Lx^ 0) plus perturbation Up = 6{x — 
xi) exp(-(x - (exp(-(y + Ly7r/2)^) + exp(-(^ - Ly7r/2)^)) , 
xi = —2Lx^ for various values of t. 



We can give some numerical evidence for the validity of the interpretation of the 
peaks in Fig. [7] as lumps in an asymptotic sense. We can identify numerically a 
certain peak, i.e., obtain the value and the location of its maximum. With these 
parameters one can study the difference between the KP solution and a lump with 
these parameters to see how well the lump fits the peak. This is illustrated for the 
two peaks, which formed first and which have therefore traveled the largest distance 
in Fig. [S] Obviously one cannot expect a true lump in this case since there are still 
some remains of the line soliton to be seen. A multi-lump solution might be a better 
fit, but here we mainly want to illustrate the concept which obviously cannot fully 
apply at the studied small times. Nonetheless Fig. [S] illustrates convincingly that 
the observed peaks will asymptotically develop into lumps. 

The precise pattern of lumps forming obviously depends on the perturbation. In 
Fig. |9] one can see the KP I solution for the same situation as in Fig. [7| except for 
the different sign of the perturbation Up ( [46| ). In this case two pairs of lump form 
first to be followed at later times by a collection of lumps. 

Perturbing the line soliton in the form uo{x) = 12sech^(x + A cos{2y /Ly))^ one 
finds the solution shown in Fig. 10 Two lumps form at the deformed line soliton 



farthest to the right. Further lumps then form at later times where the first detached 
from the deformed soliton. It appears that the perturbation develops in this case 
into a chain of lumps. The computation is carried out with = Ny = 2^^ modes 
and 10^ time steps for Lx = 16 and Ly = S. Due to the propagation of lumps of large 
magnitude for an extended time, the relative mass conservation is comparatively 
low in this case, A ~ 3 * 10~^. 
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Figure 8. Difference of the solution to the KP I equation in Fig.[7| 
for t = 4 and two lump solutions fitted at the peaks farthest to the 
right. Only very small peaks remain of these 'lumps' indicating 
that they develop asymptotically into true lumps. 



t=0.0 t=1.4 




Figure 9. Solution to the KP I equation for initial data given by 
the KdV soliton Usoi{x + 2Lx^0) plus perturbation Up = —6{x — 
xi) exp(-(x - (exp(-(^ + Ly7r/2)'^) + exp{-{y - L^7r/2)2)) , 
xi = —2Lx^ for various values of t. 
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t=6.0 



t=8.0 
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Figure 10. Solution to the KP I equation for initial data uo{x) 
12secli^(x + 0.4cos(2y/L^)) for various values of t. 



The above numerical experiments indicate that the line soliton is unstable against 
general perturbations. Whereas this might be true on long time scales, this is 
not necessarily the case on intermediate time scales. To illustrate this 'meta- 



stability' we consider in Fig. 11a perturbation of the line soliton as before uo{x^ y) = 
12sech^(x + 2Lx) + Up^ but this time with xi = xo/2, i.e., perturbation and soliton 
are well separated. The figure shows the difference between KP I solution and line 
soliton. It can be seen that the soliton is essentially stable on the shown time scales. 
The perturbation leads to algebraic tails towards positive x- values and to dispersive 
oscillations as studied in [68 . Due to the imposed periodicity both of these cannot 
escape the computational domain and appear on the respective other side. The 
important point is, however, that though the oscillations of comparatively large 
amplitude hit the line soliton quickly after the initial time, its shape is more or less 
unaffected till t = 6. 

The picture changes, however, if the code runs for a longer time as can be seen in 
Fig. [12] The soliton develops for times greater than 6 lumps after having stayed 
close to its original shape before despite strong perturbations. The solution for 
times greater than 6 remains the same within the given precision if the resolution 
is changed. For instance if we compare the result with Nt = 6400 and Nt = 12800, 
the maximal difference between the numerical solutions if of the order 10~^. The 
quantity A in the latter case is of the order of 10~^, whereas the Fourier coefficients 
decrease to 10~^. Doubling the resolution in x leads to a difference of the order of 
10~^ between the solutions. Thus one can conclude that the above solution is in 
fact correct within numerical precision. 
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t=0.0 t=0.15 




Figure 11. Difference of the solution to the KP I 
equation for initial data given by the KdV soliton 
Usoi{x + 2Lx^0) plus perturbation Up = 6{x — xi)exp( — (x — 

xi)2)(exp(-(y + L^V2)')+exp(-(y-L^V2)')), = -^cc 
and the KdV soliton for various values of t. 



Of course the periodic setting has an important effect here. Changing Ly from 8 to 
10 whilst keeping all other parameters unchanged does not affect the solution within 
the given precision. However the picture changes if we pass from = 8 to L^^ = 10 
as can be seen in Fig. [13] At time t = S the soliton is still mainly unperturbed, and 
only then the lump formation starts. This implies that the instability of the KdV 
soliton against lumps is only triggered when the former is close to the computational 
boundary. Consequently it appears that some form of 'meta-stability' of the KdV 
soliton is numerically established. For large Lx a perturbed KdV soliton can exist 
for long times, and the same could be true on R x T. But we cannot make any 
predictions on the time scales where this would be the case. 

Haragus and Pego [42 showed that the Zaitsev solution is an attractor for periodic 
(in y) perturbations of the line soliton. As we will illustrate below, the Zaitsev 
solution itself is unstable against periodic perturbations. To check the results of 
j42] one has to consider small periodic perturbations which travel with the wave. 
Above we do not have this restriction and find typically the formation of lumps that 
travel at higher speeds than the perturbed solution, thus violating the traveling 
wave condition. To get closer to the situation in [42 , we consider initial data for 
the exact KdV soliton and a small perturbation given by the Zaitsev solution uzait 
(25 ) of the form uo{x) = Usoi{x-\- L^^ 0) -\-0.1uzait{x-\- L^^ 0) with a = 1, /3 = 0.5, 
Lx = 10, Ly = 5/(5. The relative mass conservation A is of the order 10~^ and 
smaller until the lumps appear, and decreases then to 10~^. As can be seen in 
Fig. [Ml the perturbation and the soliton essentially travel at the same speed for 
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Figure 12. Solution to the KP I equation for initial data given 
by the KdV soliton Usoi{x + '^Lx^ 0) plus perturbation Up = 6{x — 
xi) exp(-(x - (exp(-(y + Ly7r/2)^) + exp(-(^ - Ly7r/2)^)) , 
xi = —Lx and the KdV soliton for various values of t. 



t=8.0 




Figure 13. Solution to the KP I equation for initial data given 
by the KdV soliton Usoi{x + ^L^^ 0) plus perturbation Up = 6{x — 
xi) exp(-(x - xi)^) (exp(-(^ + Ly7r/2)^) + exp(-(?/ - Ly7T/2)^)) , 
xi = —Lx and the KdV soliton as in Fig. 12 {Lx = 8) above and 
Lx = 10 below. 
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some time. One finds that the solution oscihates in this phase between the sohton 
and presumably the Zaitsev solution, until the instability against lump formation 
sets in. 



t=0.0 t=2.16 




Figure 14. Solution to the KP I equation for initial data given by 
the KdV soliton plus perturbation given by the Zaitsev solution, 
^0 = Usoi{x + Lj;,0) + 0.1uzait{x + L^,?/,0), a = 1, P = 0.5, 
Lx = 10, Ly = 5/(5 for various values of t. 



These oscillations of the traveling wave profile between soliton and Zaitsev solution 
is even more visible if one studies the maximal amplitude of the solution in Fig. 14 



in dependence of time as shown in Fig. 15 Three maxima can be recognized before 
the onset of lump formation. A similar behavior will be seen in the next subsection 
for a perturbation of the Zaitsev solution. In contrast to the situation shown in 
Fig. 13 the solution here does not change if a larger is chosen. The lump 
formation starts at roughly the same time as before. 



3.4. Perturbations of the Zaitsev solution. The KP I equation admits so- 
lutions that are exponentially localized in one spatial direction and periodic in 
the other. The simplest such solution (25) was found by Zaitsev jl28 l. The 
solution was generalized in |113j . We consider perturbations of t he f orm Up = 
6{x + Lx/2) ex.-p{ — {x + — y'^) which satisfy the constraint (2.1). For the 

numerical experiments we choose a = 1, /3 = 0.5, = 10, Ly = 5/(5, = 2^, 
Ny = 2^, Nt = 10000. With these settings we obtain a numerical mass conserva- 
tion A < 10~^. Notice that the unperturbed Zaitsev solution can be numerically 
propagated with the same precision as the line soliton, see [70 . 

For initial data given by the Zaitsev solution centered at —Lx/2 plus Up = 6(x + 
Lx/2) ex.p{—{x + — y^)^ we find the behavior shown in Fig. 16 One can 
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Figure 15. L'^-norm of the solution to the KP I equation of Fig. [T4j 

see that the initial perturbation develops into a lump traveling faster than the 
remaining Zaitsev solution. The latter develops in the following further lumps in a 
similar way as the perturbed KdV soliton. 



t=0.0 



t=1.4 





t=2.8 



t=4.0 





Figure 16. Solution to the KP I equation for initial data given 
by the Zaitsev solution uzait{x + Lx/2^y^0) plus a perturbation 
Up = 6{x + Lx/2) exp(— (x + — y'^) for various values of t. 
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Thus it appears that the Zaitsev solution is unstable as the line soliton against 
the formation of lumps. The precise pattern of the lumps depends on the initial 
perturbation. In Fig. 17 the same initial condition as in Fig. 16 is considered, this 
time however with a change in sign in the perturbation. The initial perturbation 
develops two lumps with further lumps forming at later times. 




Figure 17. Solution to the KP I equation for initial data given 
by the Zaitsev solution uzait{x + Lx/2^y,0) plus a perturbation 
Up = —6{x + Lx/2) exp(— (x + — y'^) for various values of t. 



Notice that the Zaitsev solution is in some sense more unstable than the line soliton. 
First the amplitudes of the considered perturbations are a factor of 10 smaller than 
the amplitudes of the solution, and still the decay to lumps is almost immediate. 
More importantly the Zaitsev solution is also unstable against a displaced pertur- 



bation {up = 6xexp(— — y'^)) as can be seen in Fig. 18 As in Fig. 11 for the 
perturbed KdV soliton, the initial perturbation develops dispersive oscillations and 
tails, but the former almost immediately destroy the Zaitsev solution and leads to 
the formation of lumps, in contrast to the situation for the KdV soliton. 

The Zaitsev solution also appears to be unstable against lump formation for pertur- 
bations with the same period in y. In Fig. 19 we show the time evolution of initial 
data from the Zaitsev solution as before, but with an amplitude multiplied by a 
factor of 1.1. This corresponds to a small perturbation of exactly the same period 
as the Zaitsev solution. It can be seen that each of the maxima of the initial data 
develops a lump-like structure. If we neglect the precise form of a multi-lump solu- 
tion and just subtract a single lump at each of the final peaks in Fig. [19) we obtain 



Fig. 20 This suggests that the formed peaks will in fact develop asymptotically 



into a periodic array of lumps. 
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t=2.8 



t=4.0 




Figure 18. Solution to the KP I equation for initial data given 
by the Zaitsev solution plus a perturbation Up = 6xex.p{—x^ — y^) 
for various values of t. 



This is however only the case if the Zaitsev solution is multiplied with a factor 
greater than 1, i.e., if the humps are enlarged as above. In this case they seem 
to evolve into lumps. If instead the Zaitsev solution is multiplied by a factor 
smaller than 1, e.g. 0.9, one gets for the time evolution of such data Fig. [2l] The 
computation is carried out for = 30, Ly = 2.5 with = 2^^, Ny = 2^ and 
Nt = 20000 time steps to yield a A ~ 3 * 10~^. In this case one again observes 
oscillations between the Zaitsev solution and presumably the KdV soliton, until 
finally lump formation sets in. 

In Fig. [22] it can be seen that this oscillatory phase can be quite long until lump 
formation starts. It appears that this oscillatory state is thus meta-stable in some 
sense, and the pattern does not change within the given precision if the code is 
run with Nt = 30000. The mechanism of the final decay to lumps is not clear. 
If one runs the code with the same settings on a larger domain {L^ = 35), lump 
formation sets in at roughly the same time. The same is true on a slightly reduced 
domain {Lx = 25). Thus in contrast to the perturbation of the KdV soliton studied 
in Fig. [l2j the proximity of the domain end in x-direction does not appear to be 
decisive. Instead it seems that numerical inaccuracies in the periodicity in y trigger 



the lump formation. In fact if the same situation as in Fig. 21 is studied with a 
smaller Ly = 1.5, the initial data show 3 maxima, which, however, develop into 
just two lumps as can be seen in Fig. [23] This indicates some symmetry breaking, 
presumably due to effects from the boundary of the computational domain in y- 
direction. The computation is carried out with = 2^^, A^^ = 2^, A^t = 20000 
with a final A ~ 2 * 10~^. It does not change within the expected limit if a smaller 
resolution in space or time is chosen. The L^-norm of the solution in Fig. 23 can 
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t=0.0 t=1.4 




Figure 19. Solution to the KP I equation for initial data given 
by the Zaitsev solution plus a perturbation with the same period 
in uo{x) = l.luzaiti^ + for various values of t. 




Figure 20. Difference of the solution to the KP I equation in 
Fig. 19 for t = 4 and five lump solutions fitted at the peaks. 
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t=0 t=1.08 




Figure 21. Solution to the KP I equation for initial data given 
by the Zaitsev solution plus a perturbation with the same period 
in uo{x) = 0.9uzait{x -\- 0) for various values of t. 
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Figure 22. L^-norm of the solution u of Fig. 



21 



in dependence of time. 
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Figure 23. Solution to the KP I equation for the situation of 



Fig. 21 with a smaher Ly = 1.5. 



be seen in Fig. 24 In comparison to Fig. 22 it can be recognized that the lump 
formation starts at a slightly later time. Thus we cannot conclude what would 
happen if the same initial data would be treated on R x T instead of T^, but it 
seems that the instabilities here could be due to a sensitivity to inaccuracies in the 
^/-periodicity. 




Figure 24. L^-norm of the solution u of Fig. 23 in dependence of time 



3.5. Numerical study of blowup in generalized KP equations. In this sub- 
section we will study the time evolution of initial data uo{x) for the generalized KP 
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equation (38). It is known [94 that the solution wih blow up for initial data with 
negative energy for p > 4/3. 

In the following we will consider the initial data 

(47) uo{x) = 60^0, exp{-a{x^ + y^)), a e R+. 

These data obviously satisfy the constraint ( |2.1[ ). 

For the supercritical p = 2 we find that these data in fact imply a negative energy 
for a = 1. The computation is carried out with = 5, = 2, = 2^^, 
Ny = 2^^ modes and Nt = 5000 time steps. The code breaks at a finite time tc 
since the solution appears to develop a singularity. It is stopped once the quantity 
A measuring numerical mass conservation becomes larger than 10~^. It can be 
checked that the Fourier coefficients at this time {tc ~ 0.0473 for the considered 
initial data) still decrease by 5 orders of magnitude. This implies that the solution 
is still accurate to better than plotting accuracy at this time. It can be seen that 



t= 0.0458 t= 0.0463 




Figure 25. Solution to the generalized KP I equation {p = 2) for 
the initial data uo{x^ y) = 6dxx exp( — (x^ + ^^)) foi" various values 
oft. 

the initial minimum develops apparently into a singularity. The precise formation 
of the latter is not isotropic in x and y as shown by the gradient of u close to the 



critical time in Fig. 26 It can be seen that the ^-derivative near the catastrophe is 
more than an order of magnitude bigger than the x-derivative. Thus it is foremost 
the ^/-derivative of the solution that diverges which then leads to a divergence of 
the solution itself. The much stronger gradients in ^/-direction are the reason why 
a higher resolution has to be chosen in y than in x. With the used parameters we 
ensure essentially the same resolution in both spatial directions. 

This behavior indicates that it is sensible to trace the time dependence of the 



norm of u and the L^-norm of Uy which are shown in Fig. 27 Both norms diverge 
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Figure 26. Gradient of the solution to the generahzed KP I equa- 
tion {p = 2) for the initial data uo{x^y) = 6dxx^^v{—x^ — u'^) foi" 
t = 0.0473. 



at the critical time. It can be seen that both norms are monotonically increasing, 
but that there is a steep divergence at the critical point. It is also obvious that the 
norm of Uy diverges much more rapidly than u. From the obtained data one cannot 
decide the exact character of the divergence near the critical time (exponential, 
algebraic). 
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Figure 27. The L^-norm of the solution to the KP I equation 
{p = 2) for the initial data uo{x^ y) = Qdxx exp(— — y"^) and the 
L^-norm of Uy normalized to 1 for t = in dependence of time. 
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For the generalized KP II equation we find for the same p and the same initial data 
that there is no indication for blowup. The time evolution of the solution can be 
seen in Fig. 
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It can be seen that the typical oscillations in x-direction due to 
the Airy term in the KP equation appear, see [68^. The algebraic fall off to infinity 
typical for KP solutions even for Schwartzian initial data is clearly visible in the 
form of tails to the left. 



t=0.0 



t=0.0343 





t=0.0686 



t=0.098 





Figure 28. Solution to the generalized KP II equation {p = 2) 
for the initial data uo{x, y) = 6dxx exp(— — y'^) for various times. 



In contrast to the KP I case, there is no indication for a blow up. In Fig. 29 the 
L^-norm of u and of the L^-norm of Uy are shown. The former seems to decrease 
since the negative hump is essentially radiated away to infinity via the tails or to 
lead to the dispersive oscillations. The L^-norm of Uy appears to approach slowly 
a finite asymptotic value. 

For the standard KP I equation {p = 1) similar initial data lead as expected to 



a regular solution as can be seen in Fig. 30 The energy for these initial data is 



positive. The solution just shows the characteristic oscillations and the tails. This 
computation and the ones below are carried out with = 2^^ and Ny = 2^ and 
Nf = 1000 time steps. The Fourier coefficients decrease to 10~^ which ensures 
spatial resolution, and the relative mass conservation is numerically satisfied to 
better than 10~^. Thus with numerical precision we can exclude the formation of 
a singularity up the considered time t = 0.15. The L^-norm of u and the L^-norm 
of Uy both decrease for large times. Thus there is no indication of the appearance 
of a singularity as expected. 



For the critical exponent p = 4/3 we obtain for the initial data (47) with a = 4 



which leads to a negative energy, the solution shown in Fig. 32 
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Figure 29. The L'^-norm of the solution to the generahzed KP II 
equation {p = 2) for the initial data uo{x^ y) = 6dxx exp(— — y'^) 
and the L^-norm of Uy normalized to 1 for t = in dependence of 
time. 




There is no indication of a divergence of the solution, and the L'^-norm of u actually 
decreases as can be seen in Fig. 33 However the L^-norm of Uy seems to increase 



without bound in accordance with the theory. It appears that no blowup happens 
for finite times, but this can of course not be decided with numerical methods alone. 
We just did not see an indication of a blowup for even longer times. 
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Figure 31. The L^-norm of the solution to the KP I equation 
{p = 1) for the initial data uo{x^ y) = 12dxx exp(— — y'^) and the 
L^-norm of Uy normalized to 1 for t = in dependence of time. 




Figure 32. Solution to the generalized KP I equation with critical 
exponent p = 4/3 for the initial data uo{x, y) = 6dxx exp(— 4(x^ + 
y'^)) for t = 0.05. 
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Figure 33. The L^-norm of the solution to the generahzed KP 
I equation with critical exponent p = 4/3 for the initial data 
uo{x^ y) = 6dxx exp(— 4(x^+?/^)) and the L^-norm of Uy normalized 
to 1 for t = in dependence of time. 
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